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Abstract 
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1 Introduction 



With a branching ratio of (1.198±0.032)% the three body decay — > e~^e~j is the second 
most important decay channel^ of the neutral pion. The dominant decay mode, vr'^ — > 77, 
with its overwhelming branching ratio of (98.798 it 0.032) %, is deeply connected to this three 
body decay. The other decay channels related to the anomalous 7r''-7-7 vertex, like e~^e~ 
and TT^ — > e"'"e~e^e~, are suppressed approximately by factors of 10~^ and 10~^, respectively. 
Another interest of the Dalitz decay lies in the fact that it provides information on the semi 
off-shell 7r''-7-7* transition form factor T^q^^* {q^) in the time-like region, and more specifically 
on its slope parameter a^^. The most recent determinations of obtained from measurements 
(2 El El of the differential decay rate of the Dalitz decay, 

= -0.11 ±0.03 ±0.08 01 

= ±0.026 ± 0.024 ± 0.0048 

= ±0.025 ±0.014 ±0.026 0, 

are endowed with large error bars, as compared to the values extracted from the extrapolation 
of data at higher energies in the space-like region, = —q^ > 0.5 GeV^, obtained by CELLO 
dl and CLEO 0, 

= ±0.0326 ± 0.0026 ± 0.0026 ^ 

= ±0.0303 ± 0.0008 ± 0.0009 ± 0.0012 

These extrapolations are however model dependent, and a direct and accurate determination 
of Ojr from the decay vr'' — > e^e~'^ would offer a complementary source of information. Let us 
mention, in this context, the proposal [HI of the PrimEx experiment at TJNAF to study the 
reaction e~7 — > e~TT^, where the neutral pion is produced in the field of a nucleus through 
virtual photons from electron scattering jlOj . Although this process concerns again virtualities 
in the space- like region, very low values of Q^, in the range well below the lowest values 
attained by the CELLO experiment, can be achieved upon selecting the events according to 
the emission angles of the produced pion and of the scattered electron jlOl ^ . 

On the theoretical side, several studies have addressed the issue of the radiative corrections 
to the decay vr'' — > e+e~7 in the past. At lowest order, the decay amplitude is of order 
0{e^). The next-to- leading radiative corrections to the total decay rate were first evaluated 
numerically by D. Joseph with the result 



"Tad ("—O 



(tt — > e~^e 7) 
r(7r0 77) 



1.0 X 10"^. 



This shows that the radiative contribution is tiny and can be neglected in the total decay 
rate. However, the differential decay rate, which provides the relevant observable for the 
determination of a^r, is sensitive to these radiative corrections. This problem was extensively 
studied in ^21 and ^Hl- In ^ cases, the two-photon exchange terms were neglected, and some 
further approximations were made (e.g. restrictions in the kinematical region and on the 
energy of the bremsstrahlung photon). Subsequently, the Dalitz decay was further discussed 

^The process tt" — > e^e~^ is currently referred to as the Dalitz decay, after R. H. Dalitz who first studied 
it more than fifty years ago [2, and who was the first to realize its connection with two-photon production in 
the emulsion events of cosmic rays. For a nice and instructive historical retrospective, see 0. 
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in connection with the omission of the two photon exchange contributions. Particularly, 
during the 1980s, the controversial question of the actual size of these contributions was 
under debate, as well as the relevance of Low's theorem in this context, cf. the articles 
quoted in ^1] . Eventually, the non interchangeability of the limits of vanishing electron mass 
and photon momentum was identified ^H] as the origin of the apparent puzzle raised by the 
contradictory results obtained previously by various authors. 

Our purpose is to provide a complete treatment of the next-to-leading radiative corrections 
to the Dalitz decay, taking into account the theoretical progresses accomplished in various 
aspects related to this issue. For instance, in the studies quoted so far, the pion was taken as 
point-like. On the other hand, the leading order amplitude involves the form factor ^^o^^* (g^) 
with virtualities up to ~ M^o ■, which is within the realm where chiral perturbation theory 
(ChPT) ^1 El E] is applicable. The details of the one-loop calculation of the 7r°-7-7* 
vertex in ChPT can be found in 19 . However, we also need to consider (among other 
contributions) the electromagnetic corrections to J^^o^^*{q'^). We are thus led to reformulate 
and extend the results described above within the unified and self-contained framework of 
ChPT with virtual photons, as it was formulated in 20^ and in [25 . Actually, it is also quite 
straightforward to include light leptons in the effective theory, as described in or, in 
the context of semileptonic decays of the light mesons, in [22]. Throughout, we shall work 
within the framework of two light quark flavours, u and d. The corresponding extension 
to virtual photons is to be found in Refs. jjl] and |25j . However, contributions of next-to- 
leading order 0{e^) to the amplitude now involve the doubly off-shell 7r^-^*-j* transition 
form factor A^o^*^*{qf,q2), but for arbitrarily large virtualities, a situation which cannot be 
dealt with within ChPT. We shall introduce and use a representation [261 127j of the form 
factor ^^0^*^* (gf , that relies on properties of both the large- A^'c limit and the short- 
distance regime of QCD. The same framework also allows us to supplement our analysis with 
estimates of the relevant low energy constants, along the lines of, for instance, Refs. |2H| and 

m- 

The material of this article is organized as follows. The general properties (kinematics, 
diagram topologies,...) of the amplitude are discussed in Section 2. Section 3 is devoted to the 
computation of the differential decay rate at next-to-leading order (NLO). Numerical results 
are presented in Section 4. A brief summary and conclusions are gathered in Section 5. For 
reasons of convenience, various technical details have been included in the form of appendices. 
Preliminary reports of the present work have appeared in Refs. |29[ I.SDj . 

2 General properties of the Dalitz decay amplitude 

In this section we describe the general structure of the amplitude for the Dalitz decay vr*^ — > 
e'^e~'-f, relevant for the discussion of the contributions both at leading order, 0{e^), and at 
next-to-leading order, 0{e^). 

2.1 Notation and kinematics 

The Dalitz decay amplitude TW^o^g+g-^ is defined as 

(e+(p+,.s+)e-(p_,s_)7(A;,A); out\7T\P); in) = i{2Ti)H^^\P-p+-p--k)M^,^,+,-^, (2.1) 

where the transition matrix element has to be evaluated in the presence of the strong and the 
electromagnetic interactions. Lorentz covariance allows to express the amplitude M.T^a^^+^-^ 
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in the form 

M^o_^^+^-^ = u{p-, s-)Ti^{p+,p^,k)v{p+, s+)£*^'{k), (2.2) 

with 

u{p-,s-)r^{p+,p_,k)v{p+,s+) = lim ie(e+(p+, s+)e"(p_, s_); out|j/,(0)|7r°(P); in) (2.3) 

given in terms of the electromagnetic current 

2_ 1 - 

j(U = -^u-fi^u - -dj^d - ipel^i^e + ■■■ (2.4) 

Invariance under parity, charge conjugation, and gauge symmetry, 

k''u{p-,s-)r^{p+,p_,k)v{p+,s+) = 0, (2.5) 

imphes a transverse structure and a decomposition in terms of four independent form factors^ 

T^'{p+,p^,k) = P{x,y)[ik-p+)p'^-{k-p-)p^]j5 

+ A+{x,y)[Ul - {k-p+hni5 - A^{x,y)[l^p''_ - {k ■ p^)j^]j5 
-iT{x,y)a^''k,^5. (2.6) 

The invariant form factors P{x,y), A±{x,y) and T{x,y) are functions of two independent 
kinematical variables, which we have chosen as (m denotes the electron mass, p^ = p\ = m^) 

{P++P-? 2 ^ ^ 1 2 

X = TTo ) < X <l, V = —pT, 



y = ^M^|^l_^)~^ ' -CTe(M» < y < ae(M», ae(.) = ^l-^. 

In the pion rest frame, these invariants can be expressed in terms of the energies of the photon 
(a;), of the positron (-E+) and of the electron (-E_) as 

1-x = 2 " 



E 

y ■■ 



E-\- — E- 



In terms of the variables x and y, charge conjugation invariance implies that the form factors 
satisfy the symmetry relations 

P{x,y) = P{x,-y), A^{x,y) = A±{x,-y), T{x,y) = T{x,-y). 

Let us note that the form factors P{x,y), A±{x,y) and T{x,y) can be projected out from 
by means of the formula 

F = Tr(A^(^_ + m)r^(|i -m)), (2.7) 



^We have omitted additional structures, proportional to A;^, which vanish upon contraction with the po- 
larization vector e*'^{k). Implicitly, wc only consider electromagnetic and strong interactions, and we assume 
that there is no P and CP violating 9 term. 
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where A^, with F = P, A± , T, are projectors satisfying k ■ Ap = 0. Exphcit expressions of 
these projectors are given in Appendix A. 

In terms of the variables x, y, the differential decay rate is given by the formula 

dr = ^^(1 - x)\M^o^,+,-^\^dxdy. (2.8) 

Expressed in terms of the form factors P, A± and T, the square of the invariant amplitude 
(summed over polarizations) reads 

polarizations 

X {[M^ox{l - y^) - Am\\P\^xMl,, - 2mP(^+ + yl_)* 

- 2mP*{A+ + A^) + 2PT* + 2P*T] 

+ 2ixM^o - 4m^)[\A+\\l + y)^ + \A^\\l - yf] 

- 8m^y'^{AXA^ + A+A*_) + 8my(l + y){AXT + A+T*) 

- 8my{l - y){A*_T + A.T*) + 8(1 - y^)|T|2}. (2.9) 

In the case m = 0, this reduces to 

lA^^o^e+e-^P = I ^'o(l - xf{M^ox{l - y'')[\P\^xM^o + 2PT* + 2P*T] 

+ 2M^ox[\A+\\l + y2) + |^„|2(i _ y2)] ^ _ y2)|^|2|^ 

As usual, higher order corrections induced by virtual photon contributions generate infrared 
singularities, even for a nonvanishing electron mass m. In order to obtain an infrared finite 
and physically observable (differential) decay rate, the emission processes of real soft photons 
have also to be considered. 

2.2 Anatomy of the Dalitz decay amplitude 

The contributions to the amplitude rather naturally separate into two main 

classes. The first one corresponds to the Feynman graphs where the electron-positron pair 
is produced by a single photon (Dalitz pair). The leading contribution, of order O(e^), to 
the decay amplitude belongs to these one-photon reducible graphs. They involve the semi- 
off-shell n^-'y-'y* vertex J^^o^^*{q^), see Fig. ^ The second class of contributions corresponds 
to the one-photon irreducible topologies. They can be further separated into the one-fermion 
reducible contributions, which represent the radiative corrections to the vr'' — > e^e~ process 
(see Fig. I2J, and the remaining one-particle irreducible graphs (Fig. OJ, starting with the two 
photon exchange box diagram, see the second graph on Fig. ^ Both types of these one- 
photon irreducible contributions to the amplitude involve the doubly off-shell 7r''-7*-7* vertex 
^^0^.^. (gf , g|). They are suppressed with respect to the lowest order one-photon reducible 
contribution, starting at the order O(e^) with the contributions depicted on Fig. ^ Let us 
now discuss consecutively these different topologies in greater detail. 

2.2.1 The one-photon reducible contributions 

The one-photon reducible topologies are shown on Fig. ^ They contain the leading order 
contribution to 7^,^0^6+6-71 and involve only low virtualities of the semi off-shell form factor 
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Figure 1: One-photon reducible diagrams 



J^^o^^*{q^). The contribution at leading, 0{e^), but also at next-to-leading, O(e^), orders 
can thus be fully treated within the framework of ChPT, extended to virtual photons. The 
general expression for this one-photon reducible part of the Dalitz decay amplitude has the 
form^ 

where 

In this and the following expressions, q = + p-. The form factor J^^o^^*{q'^) is related to 
the doubly off-shell form factor ^7^0^.^* (g^, defined as^ 

J d'xe''^{0\T{f{x)fm7r\P)) = -ie^'-"Pl^PpA^.^,^,{l\{P-lf), (2.10) 

by 

^nO-yy*iQ ) = A„Oy*y*{0,q ). 

Here the matrix element on the left hand side can be obtained by means of the LSZ formula 
from the three point Green's function (VVA), calculated within QCD+QED (i.e. with the 
QED corrections included). Furthermore, Dj^^^^q) is the transverse part of the photon propa- 
gator (the longitudinal, gauge dependent, part of the photon propagator does not contribute), 



iD^ (q) = -i - g^'^-Zg 

q^[l + U{q^)] ' 

where n(g^) is the renormalized vacuum polarization function (in the on-shell renormalization 
scheme with 11(0) = 0), and AP{qi,q2) stands for the off-shell one-particle irreducible e"'"-e~-7 
vertex function. For on-shell momenta, qi = p-, q2 = —p+, it can be decomposed in terms of 
the Dirac and Pauli form factors Fi^2{q^)i 

uK^{p-, -p+)v = u\Fi{q^)-f^ + -1-F2{q^)\a^'''qu]v, 

2m 

with -Fi(O) = 1 and -F2(0) = Og, where Oe is the anomalous magnetic moment of the electron. 
Note that the one-photon reducible part Tl['^{p^,p^, k) is gauge invariant by itself, 

A;^ri'^^(p+,p_,fc) = 0, 



Henceforth, we simply write u instead of u{p-,s-), and v instead of v{p+,s+), whenever no confusion 

3S. 
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and therefore it can be expressed in terms of form factors P, A±, T. Using e.g. the formu- 
lae (EH), (HU, (|X2|) . and (TOl) . one obtains 

^^«(...) = e'^„.,,.(.M;.,) ^^^,^^^;-^^^^,^^, iF.(xM.-.), (2.11) 

2.2.2 One-fermion reducible and one-particle irreducible contributions 

The one-fermion reducible and the one-particle irreducible topologies, shown on Fig. [21 and 
on Fig. |31 respectively, both start at order O(e^). 




Figure 2: One-fermion reducible diagrams 

Since the one-photon reducible part ■M^^^+^-^ of the invariant amplitude is transverse 
by itself, the one-fermion reducible and one particle irreducible contributions -A^^of^g+g-^ + 
■^IS^e+e-'Y together also represent a transverse subset. However, these two types of contri- 
butions are not transverse when taken separately. 

Let us first concentrate on the one-fermion reducible topology. These contributions can 
be expressed in the form 

where 

iTlf^{p+,p.,k) = (-ie)A^(p_,p_ + k)iS{p- + /c)ir^Og-g+(p_ + k,p+) 

+ iT^Oe-e+{P-,P+ + k)iS{-p+ - k){-\e)K^{-p+ - k,-p+). (2.12) 




Figure 3: One-particle irreducible diagrams 
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In this formula, 



iSiq) 



— m — T,{q) 
is the full fermion propagator, with 



(2.13) 



the fermion self-energy, while A^((7i,g2) = 7/i + and r^Og-e+(^i, ^2) are the (off-shell) 

one-particle irreducible e^-e~-7 and iT^-e~-e'^ vertices, respectively. As we have already 
mentioned, ir^'^^(p+,p_, /c) is not transverse, since A.^{qi,q2) satisfies the Ward-Takahashi 
identity 

{qi - q2) ■ Hquq2) = S'^qi) - S'\q2). (2.14) 



The solution of this identity reads 



with 



Am(9i,Q2) = Au{qi,q2) + Au(gi,Q2) 



(gi -g2) • A^(gi,g2) = 0, 



(2.15) 



(2.16) 



and the longitudinal part, which consists of any particular solution of Eq. (|2.14j) . may conve- 
niently be chosen |3J as 



+7^ 



1 - \^v{al) - \^v{^ + + ^2), (2.17) 
2 2 J 91-92 



The transverse part Aj^{qi,q2) is then parameterized in terms of eight form factors Ai, i = 
1, . . . , 8, corresponding to the eight available independent transverse tensor structures TJ^^, i = 
1, . . . , 8 (we will not reproduce them here, for a detailed discussion and explicit expressions, 
see Refs. |SI] and |S2]) 



A^(9i, 92) = Yl ^^('^1' il - 92)')r;(9i, 92) 



(2.18) 



The decomposition ()2.15|) of the vertex function A(9i,92) induces the corresponding decom- 
position of r^'^^(p+,p_, /c). 



(2.19) 





WW 



Figure 4: One-photon irreducible contributions at one-loop level. 
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with 

k^^T]f''-'^{p+,p^,k) = (2.20) 

and 

k^'rlf^'^{p+,p_,k) = er^Oe-e+{p- + k,p+)-eT^Oe-e+{p~,P+ + k) 

-S'^{p^)S{p- + k)eT^o^-^+{p_ + k,p+) 

+eT^o,-MP-,P+ + k)S{-p+ - k)S-\-p+). (2.21) 

Therefore 

uk ■ T^^^'^{p+,p.,k)v = eu[Oe-e+(P- + k,p+) - r^Oe-e+{p^,P+ + k)]v. (2.22) 

This non transverse piece should be cancelled by the contribution u{k ■ T^^^)v of the one- 
particle irreducible graphs. In addition, the transverse part T]f^''^{p^,p^,k) admits a rep- 
resentation of the type 1)2. 6|) . with appropriate form factors F^^^''^ {x,y), A^^''^ {x,y), and 
T^^^''^ {x,y), up to possible terms proportional to k^, which cancel when contracted with 
e''{k)*. 

As for the vertex T^0e-e+{Q2, Qi), it can be decomposed (using Lorentz invariance, the 
Dirac structure of the inverse fermion propagator S~^{q), and charge conjugation invariance) 
as 

^■K0e-e+{Q2,qi) = PnOe-e+iQhQih^ + l5S^'^{-qi)A^0e~^+{q2,ql) 

+ S-\q2)l5A^0e'e+{qlql) + S-\q2)l^S-\-qi)T^,,-,+ {qlql), (2.23) 

where Pj^Oe-e+i ^7r0e-e+; ^'^^ '^■K°e-e+ 3-^6 Scalar form factors, which, as a consequence of 
charge conjugation invariance, satisfy the additional relations 

The form factor ^,^06-6+ is then related to the on-shell vr^ — > e+e" amplitude, 

M^O_e+e- = uiP-,S-)-/5v{p+, S+)P^Oe-e+ (m^, m^). 



In terms of the form factors (I2.2.SI) we can write 



= ee^{kru{A^{p^,p^ + k)Sip^+k)P^o,-,+ {m^ + 2{k.p^),m^h5 
+PnOe-e+{m^ + '^{k ■ P+)h5S{-p+ - k))A^{-p+ - k,-p+) 
+K^{p-,p- + k)-i^A^Oe-e+{m^,'rn^ + 2(A: -p-)) 

+75A/.(-p+ - k, -p+)A^o^-^+{m^,m^ + 2{k • p+))}v. (2.25) 
At leading order in the fine structure constant a, F^Og+g- is given by 

{P + iO) [(gi + g2 - /)2 + zO] ^2 - / - m + iO ' ^ ^ ^ 
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where -4.^o^*^* (gf , is now restricted to its pure QCD part. The corresponding expression 
of Tjf^{p^,p^, k) then reads 



(Z2 + iO)[(P-02+iO] 

_ + |^-m + iO ( -m + iO 



1 1 

'^1' '2 71. ^ I \n 7. Til. 1 ™ I ;n ^'^ 



:^ - / - m + iO ^'^ + ^ + m + iO 



(2.27) 



The general properties of the form factor A^o^^ (qf, gl) are summarized in Appendix B. Here 



we only note that the short distance behaviour of A,^o^*^*{P , (P — /)^) in QCD makes it act 



as an ultraviolet cut-off, so that the loop integral on the right-hand side of (|2.26|) actually 
converges. 

Finally, the one-particle irreducible part of the amplitude 

Ml?i^^^.^ = urjf'{p+,p.,k)ve^{kr 
starts at the order with the box diagram of Fig. 



ir]f'{p+,p^,k) = eV™/^|^Aor7*(^''(^-0') 



{P + iO)[iP-iy + iO] 



^^p^ — r — ~^^^7r~j — ) z^^-' (^-^^^ 

— I — m + i\) p_ + f — I— m + i\) 

which is also ultraviolet finite. At this order, one verifies that the sum uT]f^{pj^,p-, k)v + 
urjf^{p^,P-,k)v is indeed transverse. 

3 The NLO differential decay rate 

The leading order amplitude corresponds to the 0{e'^) one-photon reducible contribution, 
evaluated at lowest order in the extended chiral expansion, i.e. with = 1, F2(g2) = 

11(5^) = 0, and with the form factor ^^o^*^* (/^, (P — /)^) reduced to its expression for a 
pointlike pion, i.e. a constant, ^^o^,^, = — Nc / I^tt"^ F^^ , fixed by the chiral anomaly. The 
leading order expressions of the form factors P, A± and T are then given, for Nc = 3 and 
according to Eq. by 

= 0, 



47r2F^M2 



X 



Note that in the limit m ^ 0, only the form factors A^ survive. The square of the leading 
invariant amplitude summed over polarizations is, according to (|2.9j) . 

6 n _ 



l-^^o«-.e.e-/ = [^.^ox(l + y^) + W] (3.2) 
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and the corresponding partial decay rates read 

dr^o _ M^o (1 - 
dxdy ~ (47r)4 

dr^o ^3 gM^^ (1-x) 



[M^oxil + y^) + 4m% 

3 

- aeixM^o) (xM^o + 2rr?). (3.3) 



dx (47r)4 3 

The next-to-leading corrections to the differential decay rates will be described as 

dr dr^*^ 



dxdy ' dxdy' 

dr , dv^o ^ • > 

d^ = '^^")^- 

Knowledge of the corrections 5{x, y) and 5{x) to the Dalitz plot distributions allows to extract 

information on the QCD part of the form factor T^o^^* {q^) from the experimentally measured 
decay distribution. For instance, if the form factor is approximated by a constant plus linear 
term 

•^ttOtt* (9^) = ^-Ko^r (0) [1 + +■■■], (3-5) 

the slope parameter a,r is obtained from 

^ - Sqed{x) ^ = ^ [1 + 2^ ^'t], (3.6) 

where the QED part, 5qe_d(x), of the corrections S{x) will be specified below. 

For the purpose of the following subsections, wc introduce functions 6'^{x,y) and (5*(x), 
i = IjR, IjIR, ■ ■ ., measuring the magnitude of various 0{e^) and/or 0{e^p^) corrections dF' 
to the leading order decay rate. In terms of the corresponding corrections to the invariant 
amplitudes, T{x,y) = T^'^{x,y) + {S''T){x,y), etc., one has^ 

6^{x,y) = -4^' n^2,^, , xlm{8m(<5-r)(x,y) (3.7) 

+ [M^ox{l + yf - 4m2](5M+)(x, y) + [Ml,x{l - yf - 4m2](<5M_)(x, y)} 



and 



SUx) = - — z / dy— — ^^5 S'(x,y). (3.8) 



3.1 NLO one-photon reducible corrections 

The computation of the corrections belonging to the one-photon reducible type of topology 
requires the evaluation of several quantities beyond leading order. Thus, next-to-leading cor- 
rections, of orders 0(p^e^) and 0(p^e^), to J^^Q^^*{q^), as well as corrections of orders 0{p^e^) 
and O^p^e^) to the electromagnetic form factors Fi{xM'^q), F2{xM^q), and to the vacuum 

®As usual, the NLO corrections to the decay rate arise from the interference between the leading and NLO 
amplitudes. This explains why there is no contribution involving {5^P){x,y) in these expressions, given that 
P^'~'{x,y) vanishes. 
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polarization function Il{xM^o), have to be evaluated within the framework of (extended) 
ChPT. These corrections involve one loop graphs with virtual pions, photons and electrons, 
and local contributions given in terms of counterterms. The interested reader may find the 
details of these calculations in Appendices lUl and IdI The corresponding NLO corrections to 
the Dalitz distribution read 



2Re 



+ 



and 



2Re 



2xM\ — , 9 . 



1 



a^l^^ (xM^o ) -Ti{xMl, ) + Fi (xM^o ) - 1 + 



2M^„x + 2m- 



-F2{xMl,) . (3.10) 



(3.9) 



The expressions of the various quantities appearing in these formulae are displayed in Eqs. 
(frr2T]) . (|m)) . l|TT6)) and ipTf)) of Appendix O and El Let us just mention here that at NLO 
the Dirac form factor -Fi(s) develops an infrared singularity, 



IRdiv 



m 



2ir 



1 + (s - 2m2 



1 



In 



1 - CTejs) 



+ ivr 



where is a small photon mass introduced as an infrared regulator. Thus, the infrared 
divergent part of the one-photon reducible corrections reads 



div 



In 



m 



1+1 



2m^ 



1 



1 - aejxMy, 
xMlJ ae{xMl,) "Al + (^e{xMl,] 



In 



(3.11) 



3.2 One-photon irreducible contributions 

The evaluation of the contribution 5^'^^^{x,y) involves the QCD form factor ^7^0^*^* (g^, g^) 
for arbitrary virtualities. This in turn addresses non perturbative issues beyond the low 
energy range covered by ChPT. While the asymptotic regime can be reached through the 
short distance properties of QCD and the operator product expansion [331 EU) there still 
remains the intermediate energy region, populated by resonances at the 1 GeV scale, to be 
accounted for. If one restricts oneself to approaches with a clear theoretical link to QCD, 
the large- A'^c' framework is almost the only available possibility^. In Refs. |27| I26j. the form 
factor ^^0^.^. ((7^, gl) Las been investigated within a well defined approximation to the large- 
Nc limit of QCD, which consists in retaining only a finite number of resonances in each 
channel. Details of this approach, as far as the form factor ^^^o^*^* (g^, gl) is concerned, are 
to be found in Appendix El Thus, upon inserting the expression of Eq. ()B.4|) in the form 



ALMD{l\{P-lf 



3M^ 



Pf 



M2 



P[{l-Pf 



Ml 



+ 



2M^ + KV 



(P-Mmi-py 



M2] 



^Large scale numerical simulations on a discretized space-time might become an alternative in the future. 
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into ()2.26|) and l|2.28() . the integral over the loop momentum can be done and expressed in 
terms of the standard one-loop functions Bq, Cq and Dq defined in AppendixEl It is, however, 
much easier, and equivalent^, to proceed within the framework of an effective lagrangian 
approach using ChPT with explicit photons and leptons j23j . that we now briefly describe. 
Thus, we take in (|2.26|) and (|2.28|) the leading order constant expression 

A^?^,^, = -NJUtt^F^. (3.12) 

This is a good approximation only for sufficiently low loop momentum, P <^ A'jj , where Ah ~ 
1 GeV is the hadronic scale typical for the non Goldstone resonance states. The intermediate 
and asymptotic ranges are, however, not treated properly using this effective vertex. One of 
the consequences is that the loop integral 1)2.26(1 with this constant form factor is ultraviolet 
divergent. Within the framework of an effective low energy theory, the difference between 
the exact and the low energy effective vertex can be taken into account by a counterterm 
contribution stemming from the Lagrangian [20] 

Cpe-e+ = § [^yi^i^7'ij[xi{Q\D^UU+-D^U+U) + X2{U^QD^UQ-UQD^U+Q)]. 

Let us note that these counterterms are also necessary to cure the ultraviolet divergence that 
arise in the loop integral of Eq. (|2.26|) with a constant form factor. /2pe-e+ generates a local 
TT^{p) — > e~{q)e~^{q') vertex of the form 

or, in terms of the decomposition l|2.23() . 



pCT 




. 8mx 
= —I 


aCT 




Ax ( 


rpCT 




= 0. 



(3.13) 

In the above formulae, x stands for the relevant combination of the effective couplings,^ 

(3.14) 



x = -^^^ = x'if^) + 3 



11/, 

-(In 47r - 7 + 1) 



d-4 2' 



which decomposes into a finite, but scale dependent, renormalized part x^il^)^ ^-iid a well 
defined divergent part [H^. We therefore split the amplitude -A/i^ol^e+e-^ parts 

7r^ — ye^e 7 tt'-* — >e^e 7 tt'-' — >e^e 7' 

corresponding to the loop (computed with the constant form factor ^^o^*^*) and counterterm 

contributions, respectively. Because -M^o^g.^-^ is gauge invariant, it can be decomposed into 

^The results obtained within the two approaches will differ by terms of the order 0{'m? /My). 
*We use the convention where both the loop functions and bare couplings are renormalization scale depen- 
dent - see also Appendix El 
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form factors according to (|2.6() . with 



,8mx ( o. 



\ttJ M4o(l-x)2(l-y2)^ 



_2mx /ct\2 



WJ M2,(l-x)(l-y2)' 



For the corresponding decomposition 6^'^^^{x, y) = 5^'^^^' '^^(x, y)+5^"'^^'' y), one finds, 

upon using Eq. (|3.7j) . 



rin?x 



7r7 (1 - x) (1 - y2) [m2ox(1 + y^) + 4^2] 



(3.15) 



The interference term of the loop amplitude with the lowest order one-photon reducible am- 
plitude A^^o^_^g+g-^ results in 



a 



.-kJ [x(1 + y2) + z,2] (i_^)2 
xRe { i(x - lf{{x - 1)2(/ - 1) - 4z.2y2)M4,Z)o 

(X - 1)((X -l)(y- l)(2/2 + _ X - y - 1) - 4zy2(y2 - y + 1)) 



+ 



+ 



4(y-l) -^'^-o 
(x - l)((x -l){y + l){y^ + l){xy + x-y + l)- Au^ {y^ + 2/ + 1)) ,^2 

4(y + l) 

(x - l)((z/ - l)2((x - l)2(y2 + 1) + ^2(^ _ 2)) + ^4) 

4(y - 1) 

(X - l)((y + l)2((x - 1)2(1/2 ^ 1) ^ ^2(^ _ 2)) + ^4) 

4(y + 1) 



; ((x - l)(y - l)2(x(3y - 5) + 2) + v\y - l){x{y - 4) + 3) - j.^) 



8(y-l)2 



ml 

, ((x -l){y + l)\x{?,y + 5) - 2) - t/2(y + l)(x(y + 4) - 3) + u^) 

8(y + l)2 ml 

2(x-l)2(y2_ 1)2(3.(7 + 3y2)_iQ) 

v'^{x - l)(y^ - 1)((4 + x)y2 llx - 16) 
-u\l+y\x{l + y'')-'i)+v\l+y'')) 



Bo (3.16) 



16(-1 -Fy2)2 



rU 

2 2"-°0 

mZ.m 



-|- (x — l)x -|- 



5z^2(a;(3 + y2) _4) j,4((g^_8)y2_^7^ 



l/^X 



2(y2_i) 4(x-l)(y2-l)2 



32(x - 1) V (y + i)"i+ (y - i)"i- 



:)}• 
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In this formula we have used the shorthand notation 

^0° = 5o(0,m2,m2), = Bo{ml,0,m^), 

Do = Do(?n.^O,m^M^o,m^ ,m^,0,m^m^O), 

where 

ml = m^ + -{l-x){l±y)M^o, 

and Bq, Cq and Dq are the standard scalar loop functions (bubble, triangle and box) listed 
in Appendix El Both 5^"'^^'^'^ {x,y) and 5^"'^^''^°°^{x,y) contain divergences, in the form of 
poles at d = 4, contained either in the bare counterterm or in the loop function Bq. Prom 
Eqs. (|E.2|) and (|E.3|) . one deduces 

= -48 -i- (^) Jfj (3.17) 

a — 4 Vvr/ (1 — x) (1 — y^) [M^qX( \ + y^) + 4771"^] 

As follows from Eqs. 1)3.14(1 and ()3.15|1 . these divergences cancel in the sum, so that 5^'^^^{x, y) 
is both finite and independent of the renormalization scale /u. 

In the literature, the explicit calculation of Ai^Z^^ . _ = Ai^t^ + _ + M.'^o^ + _ , 
when considered at all, was discussed in the approximation m = and assuming the pion to 
be pointlike, cf. [H], i.e. ^^0^6+6-7(^^1 {P - 0^) = -^^o^e+e-7' see (jXT^ . Let us note that 
in this case the ultraviolet divergent part of A^^ol^e+e-7 ^^'^^^'^^s- Indeed, the divergent part, 
for m = 0, is contained in the expression 

„ / 1 



2ik.p_) ) + ( 2ik"p,) )hP-Kiky (3.18) 

Upon using the identity 

e^""^/aP/37^/7" = 2i[/ (/ • P) - ffh^ 

and an effective substitution l^lp — > CPgaf3 (where C depends on the cut-off prescription 
used to regularize the divergent integral), one obtains 



^^5'-fe+e-7;.» = -2e^^^o^_e.e-,(C - 1) I 



dH / 1 



2[k-p^) ) V 2(A:-p+) 
The two terms in the curly brackets cancel each other as a consequence of the identities 



u{p-)f{i>^ + ]^) = 2{k-p+)u{p^ 
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and thus the ultraviolet divergences are absent in the limit m — > 0. This limit appears at 
first sight to be a very good approximation, because the relevant dimensionless parameter, 
= (2m/M^o)^ ~ 5.7 x 10~^, is tiny. This indeed turns out to be the case as far as the 
corrections to the total decay rate are concerned. However, when considering the differential 
decay rate, this simple argument can sometimes be misleading, as we discuss in the following 
subsection. 



3.3 The Low approximation 

Let us first briefly comment on the possible approximation of the above result by the appli- 
cation of the Low theorem. Since we are dealing with a radiative three body decay, we can 
borrow from general results [HE] and obtain 

M^0^,+,-^ = M^SZe+e-y + ' k)), (3-20) 

with 

A<iC...-, - ee,W-.[ ^^--"<^-^-^;|^-:g'-'-"^"""°'' (3.21) 

where the important point is the absence of contributions that are independent of in the 
difference - ■^^oZe+e-'y (^^^ Appendix|Fj. 

Note that the on-shell amplitude P^Og+g- (m^, m^) , evaluated to the order under consid- 
eration, reads 

P^0g+g-K,m2) = +PfoI+g-K,m2), 

where 

/^7+^4"^^"^') = -i(^)'£[5 + 3i?o(0,m^m2)-M2oCo(m^M^2o,m^m^O,0)]. 

This means, using (|F.2I) and (|3.7I) . that the Low amplitude ■M^oZe+e--y geiisrates the following 
correction 

, , , ^ , i^'^x 1 

{x,y) 



X Re[5 + 2x + 3Bo{0,m\m^)-M^oCo{m\M^o,m\m^,0,0)], 

which corresponds exactly to the single pole part of the complete one-photon reducible am- 
plitude M^l^^+ _ = M^V^''+°^ for X ^ 1. When integrated, the Low 
contribution to 6{x) becomes 

nJ {2x + u^){l-x)ae{xM^,) \l - ae{xM^,) 

X Re[5 + 2x + 3Bo{0,m^,m'^)-M^oCo{m^,M^o,m^,m^,0, 
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Notice that S^°^ is suppressed by the factor u'^ and vanishes in the hmit m — > 0. This was 
the argument for the conjecture that in this hmit 8^"'^^ does not develop a pole when x ^ 1 
|14j and the contributions of I'yIR topologies can be safely omitted. In fact this conjecture 
is not quite true for several reasons we shall briefly discuss now. 

First, the Low correction is not numerically relevant for almost the whole phase space. 
Because of the suppression factor z^^, the corrections 5^"^{x, y) and 6^°^{x) become important 
only in the experimentally irrelevant region where 1 — x ~ (when y is fixed), or where 
1 _ y2 ^ (^-^j^gji X ^> i^^ is fixed), i.e. for \y\ ~ yma.x{x) = y^l — v"^ /x. This is in fact 
no surprise, because it is precisely this corner of the phase space where Low's theorem is 
applicable. Indeed, the standard textbook derivation of the Low amplitude involves (and 
assumes the existence of) the power expansion of the form factors corresponding to the off- 
shell T^^e'^ {qi)e~ {q2) vertices r^Oe-g+ (^2) 9i) and rjrOg-g+ ((72) 5i) in powers of at the points 

= rn?, where qi = qi + k. This means, that the relevant expansion parameter is 

A± = ^ - 1 = — 2— = ^ 1 - 2; (1 ± y)- 

Therefore, the 0{k) terms in the formula 

M^0^,+ ,-^ = M^SZe+e-, + 0{k) (3.22) 

are small for A± <^ 1, and not just for 1 — x <C 1, as one could naively expect. 

There is another subtlety connected with such an expansion. According to Low's theorem, 
in the region of its applicability one would gather from 1)3. 22() 

5i^^^(x,y)-5^«-(x,y) = 0(l) 

with the 0(1) term independent of k (recall that the leading order amplitude -A^^o^g+g-^ is 
of the order 0{k)). However, the points = m? do not belong to the domain of analyticity 
of our e'^ {qi)e~ {q2) amplitude, because of the branch cuts starting at qf = m? due to the 
intermediate 6^7 states. As a result, the asymptotics of the amplitude for x ^ 1 will also 
contain, apart of the pole terms, non analytical pieces, like non integer powers and logarithms. 
This means we can expect 

7W^o^g+g-. = M^SZe+e-y + 0{k\nk) + 0{k) + ... 



rather than (|3.22|) and, result 

^^^■^^(x, y) - 6^"'"{x, y) = 0(ln(l - x)) + 0(1) + . . . . 

The Low correction therefore does not saturate the singular part of 6^'^^^{x,y) for x ^ 1. 
This can be verified explicitly at lowest order. Using the asymptotic form of the loop functions 
(cf. AppendixEl we find from (|3.16|1 . for A± <^ 1, 

6'^^^ix, y) - 5^-(x, ,) = (^) [((1 - y? - ln(^(l - x){l - y)^) 

+ ((1 + yf - v^) ln(i(l - x)(l + y)^)] + 0{l) + .... 
\1 /J 
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To conclude, the Low amplitude does not provide us with a numerically relevant estimate of 
S^f^^(^x,y) in the kinematical region of interest. 

On the other hand, for A± ^ 1, which is satisfied practically in the whole relevant 
domain of xand y (with the exception of the region where a; ~ 1 — z^^ or |y| ~ yma^ix) = 
■sJ^V^^i^^^Jx with X ^ z^^), we can approximate the correction 6^"'^^{x,y) with its m ^ 
(x < 1 and \y\ < ymaxix) fixed) limit with very good accuracy. Note that in the case m = 
the loop integration is infrared finite for k ^ 0, and the ultraviolet divergences as well as the 
counterterm contributions vanish. Using the corresponding asymptotic formulae for the loop 
functions (see Appendix El i the limit can be easily calculated with the result (cf. |15j where 
this approximative formula was published for the first time) 



X, 



X 



+ 



(l-x)(l+y2) 6 



In 



1 



Li2(l--(l-x)(l-y) 



l-x)(l-y))ln(^(l-x)(l + y) 



(3.23) 



Li2(l-^(l-x)(l + y) 



Notice the presence of the ~ (1 
expression (jXTTHl for 6^^^^{x,y) 



term, which stems from the following part of the 



"pole 



ix,y) 



a 



ttJ [x{l + y^)+u^] (1 



Because the limits m ^ and x ^ 1 are not interchangeable, as pointed out in Ref. ^S], for 
m 7^ the contribution of such a term is cancelled by the expansion, in powers of (x — 1), of 
another term, namely 



a 



u^x 



ttJ [x(1 + y2) + u^] 32(x - 1) \iy + l)ml {y - l)ml 

so that the only pole terms which survive are suppressed by a factor i^^ according to Low's 
theorem. In the limit m ^ we also obtain |15j 



«'-'"'(x)|„_„ 



-) 

IT/ 



ln2(l -x) + 



2x 



1 



■ ln(l — x) 



x(2x2-3x + 3) /TT^ ^ ^ ^ 

+ (i-.P H it-'-''"' 



x(5x + 3) 



4(1 -x)2 



(3.24) 



which provides an excellent approximation to the exact m ^ result in the whole relevant 
range of x. 



3.4 Soft photon bremsstrahlung 

The virtual photon corrections described in the previous subsections produce infrared di- 
vergences, which were regularized by introducing the soft photon mass m^. As usual, an 
infrared finite result is obtained at this order upon adding to the decay rate the real photon 
bremsstrahlung correction. This corresponds to the radiative process vr'' — > e"'"e~77B, which 
cannot be distinguished from the non-radiative one for energies of the bremsstrahlung photon 
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smaller than the detector resolution AE. In the soft photon approximation, the amplitude of 
the radiative decay is related to the leading matrix element by 

u . ( P--^*BikB) p+-e%{kB) \ LO 

\ -kB p+-kB J e e 7 

where kB and £^(^5) are the momentum and polarization vector of the bremsstrahlung 
photon, respectively. Squaring the amplitude and summing over polarizations, one obtains 



ip+ ■ kB){P- ■ kB) {p+ ■ kB? {p- ■ kB)\l --"^----r 

The corresponding correction appearing in (||-{.4j) is then 



d^kfi / 2{p+-p^) -2 _2 



|kB|<A£' 



m m 



(2^)32fc0 V(P+ • kB){p- ■ kB) {p+ ■ kB)^ {p- ■ ks) 



where k^ = yk^ + The correction 5^{x,y) can be expressed, in terms of the standard 
integral 

d^kB 1 



J{q,q') 



(27r)32A;0 {q ■ kB){q' ■ ks)' 



as 



S^{x,y) = ['2{p+-P-)J{p+,P-)-ra^J{p+,P+)-rn?J{P-,P-)) 

= (? {{xM^o - 2m2)J(p+,p_) - m'^J{p+,p+) - m^J{p^,p^)) . 

Let us note that the integral J(g, q') is not Lorentz invariant and the result is therefore 
frame dependent. On the other hand, the infrared divergent part of J(g, q') is given by the 
invariant expression 



, , 1 , /4Ai?2\ /■! dx^ 



'12 



2(27r)2 "^7 y Al/2(s,g2^g'2) " L _ g2 _ _ ^l/2( 



f4AE^ \ 1 ^Js-q^-q'^ + \^^^{ s.q\q'^) \ 

s,q'^,q'^) J ' 

where A(x, y, z) = x"^ + y"^ + — 2xy — 2xz — 2yz is the triangle function and s = {q + q')"^ ■ 
The infrared finite part can be transformed to the form 

, ^ r V + (1 - 2;)^'° 

JiR fin[q, q ) 



2(27r)2 7o [xq + (1 - x)q']'^ \xq + (1 - x)q'| 
xq^ + (1 — x)q"^ + |xq + (1 — x)q'\ 



X In 

^xq^ + (1 — x)q"^ — |xq + (1 — x)q'\J 
In an arbitrary frame we can easily obtain 



2(27r)2 g2 



mf^U^in^'^^-l^l 



\ m'l J |q| + 
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For q 7^ g', the calculation of the explicit form of J{q, q') is much more complicated. In the 
center of mass of q and q' , with q^ = q''^ = m?, the integral Jjr fin{Q, q') simplifies considerably 
and we obtain 



1 



+ - In ' ^ ^ 



+Li2 
J{q',q 



l + cj(s) 



Lio 



2(27r)2 m2 



In 



2 



as \l-as) 



where a{s) = (1 — Am'^/sy^'^ and X2{x) = ^[^^2{x) — Li2(— x)] is the Legendre chi-function. 

If we interpret AE as the photon energy resolution in the center of mass of the Dalitz 
pair, we find 



with 

CiR div 

and 



/ 4AE^ \ 



2m^\ 1 n+aejxMX 



(3.25) 



a 



1 



2m^\ [1^ [l-ae{xMl,] 
- m ' 



l + (Te{xMl,] 



In 



xM^ 



+Ll2 ( ^ ) - Ll2 ( ^ ) - 4x2(fTe(2;M^o)) 



+ In 



l + ae{xMl,) 
l-ae{xMl,) 



Summing divi^iV) ^N'LO^^^yhR div, as given by (|H.25j) and (|3.11|) . we explicitly 
achieve the expected cancellation of the infrared divergences, 



^Fr div{x,y) + S]^l^o(^,y)iR div = ^ 



2m^ \ 1 ^/ l-ae{xMy, 
xMlJ UeixMl,,) V 1 + ae {xMl„ ] 



4 Numerical results 



In the previous sections we have classified the NLO corrections according to the general topol- 
ogy of the corresponding Feynman diagrams. The complete correction SNLo{x,y) (see (j23I) 
is then given by the sum of the individual contributions of the one-photon reducible, brem- 
sstrahlung and one-photon irreducible graphs: 

SNLo{x,y) = 4fo(x,y) +5^(x,y) +5^^^«(x,y). 

A similar decomposition holds for 6nlo{x)- The resulting formulae contain several renormal- 
ization scale independent combinations of the a priori unknown low energy couplings and of 
chiral logarithms. 
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4.1 Inputs 

The contributions of the low energy couphngs to 5NLo{x-,y) and Snlo{x) are contained in 
^i7/R;CT^^^ y) and in a%\^J{xM^o), see for instance Eqs. (PTT^ . and (fTTnl) . We define 
the scale independent quantities 

2 1 A^2_ 



\ 77-± / 



+4x - 



TtJ (l-x)(l-y2) + y2) + ^2] 

2 1 /\/f2 



+ ^UJ (l-x)(2x + Z.2)^^(^M2o) ''Vl-<^e(xM2„);' 

which contain the contributions of the low-energy constants, with 

x = x (^) + 2!"-^- 

The differences ^'^fx)^ = ^nlo — ^mx) ^'^^ expressed in terms of the known physical observ- 
ables (and the detector resolution /S.E) and represent therefore a numerically unambiguous 
part of our calculations. 

For the combinations Ci, C2, and ICp, we have reasonable estimates based on resonance 
approximations, as described in AppendixlO For x the LMD approximation was studied in 
|2n] , with the result 

11 . .Fj 3, /M2^ 



xW/^) = ^-4vr^--ln^^, 
For numerical calculation we take the LMD values, with My = 770 MeV: 

Ci = (2.2 ±0.3) X 10"^ 

Kf = -28 ±8, 

C^{n = Mv) = (1.5 ±0.5) X 10-\ 

x''{^^ = Mv) = 2.2 ±0.7. (4.1) 



4.2 Radiative corrections to the differential decay rate 

The traditional point of view is to separate from the complete NLO corrections the pure 
electromagnetic part Sqed, which includes the I'jR graphs with the virtual fermion and 
photon loops only and bremsstrahlung contribution, together with IjIR diagrams (the latter 
were usually omitted in the analysis of the experimental data we shall comment on 

the consequences of this omission below): 

Sqed = S^^%,^ioops ± ± <5i^^^, (4.2) 
where, cf. Eqs. (|(L21|) and (|D.1|) . 

6'^%,^ioops{x)=d'^'^{x)-2Re 



a%fo^{xM-^o)-n^±{xM-^o) 



(4.3) 
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Following this point of view we present here separate plots^ for (^]vxol7:^'ooP'5(^) ~^ 
where the experimental situation is parameterized by the detector resolution AE (for which 
we take AE = lOMeV, 15MeV and 30MeV, see Fig. IHl), and for 5^'^^^{x,y) together with 
^iiiRi^x) (depicted in Fig.iJ. 



.08 




0.2 0.4 0.6 0.8 1 



Figure 5: The traditional QED corrections (without I'jIR contributions) for different detector 
resolutions AE = lOMeV (solid curve), 15MeV (dashed curve) and 30MeV (dotted curve). 
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Figure 6: The one-photon irreducible corrections (triangle and box diagrams): the dashed 
line represents 6{x,0), the dotted line S{x,ymax), and the solid line on the right 6{x). 

In the latter two we use the value of x^iMv) mentioned above. As we have discussed in 
the previous section, 6^'^^^{x) can be safely approximated with its m — > limit (|3.24)) for 
almost the whole range of x; the same is true for 6^'^^^{x,y) for \y\ < yina.x{x), the difference 
between 6^'^^^{x,y) and (|3:23l) can be seen for y ~ ymax(a^) in Fig. [3 From these figures we 
can conclude, in agreement with |I3], that the usually neglected IjIR corrections 6^'^^^{x) 
are in fact important particularly in the region x > 0.6, where they are in absolute value 
larger than 1% (up to ~ 4% for x ~ 0.9), and comparable with S]^I^q\^^^ loops + the same 
is true for S^'^^^{x, y), which is almost independent on y (except for a very narrow region near 
\y\ ~ UmsLxix), cf. previous section). The complete pure electromagnetic corrections 5qed{x) 
and SQED{x,y) are represented in Fig. |SJ 

Let us now change the point of view a little bit, and split the corrections in the way we 

^We do not show S]^^q\^^^ ioop3{x,y) + S^(x,y), because the y dependence is suppressed here by the factor 
i/^ for X > , and is therefore negligible in the relevant region of x. 
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Figure 7: The difference between 5^'^^^{x, Umax) (dotted line) and 6^'^^^{x, y, 
line). 



max J I m- 



+0 (dashed 



have described in the beginning of this section, namely 

r rknown , s:LEC 

ONLO = + O^j^o- 

Here (J^nown (;jiffgpg from 5qed (with x set to zero) by the contribution of II^±{xM'^q) as well 
as by corrections stemming from chiral pion loops, 



rknown , 



QED(X,y)\x=0 



+ 2Re 



where 



ChPTn2\ 
"•NLO ) 



16^2^2 



In 



M2, 



M2„ 6M2, \3 M2 



1 M2± 

- + ln- " 



167rV^+(/^)J,+ (/^ 



(4.4) 



Separation of the numerically unambiguous part S^'^°'"'^ allows, at least in principle, to con- 
strain the relevant combinations of the chiral low energy constants C[, C2 and ICp from 
experiment. Fig. IHl shows both (J^nown ^^^^ Snlo, which allows to appreciate the effect of 
the counterterms. The difference (5'^'^°™"(x) — 6qed{x) is particularly important for x > 0.7, 
where it represents a correction larger than 1%. 
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Figure 8: On the left hand side, the complete electromagnetic corrections Sqed with the same 
detector resolutions as in Fig. [HI The right hand side shows 6QED{x,ymax) (dotted curve) 
and 6qed{x,0) (dashed curve). 
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Figure 9: The NLO correction 5^^°'"^{x), without counterterm contributions (left), and the 
complete correction SNLoi^) (right). The error bars show the uncertainties coming from the 
counterterm determinations. The detector resolutions are as in Fig. [HI 



Let us split further 

eLEC _ >:LEC , ^LEC 
(^NLO — (^LMD + "QED^ 

where Sj^f.fjj = S^^q\!Cp=o is the part for which we have theoretical prediction based on the 
LMD approximation; S^'^"^'^ + S^ffj^i is shown in the Fig. 121 the error band stems from the 
estimate of the uncertainty of the LMD values ()4.1|) . 

For completeness, we present our theoretical estimate of Sq^^, based on (|(].24|) (cf. (|4.H) ). 

SQED{^,y) = 5qWd{x) = ^}Cf = (-8 ± 2) X W'^. 

Note that this value is comparable with the estimated uncertainty of S^ff^ arising from the 
uncertainties in Ci and C2(/^ = My) as given in Eq. (|4.1|) . 

4.3 Decay rate 

Let us make a brief comment on the total decay rate. At the leading order we reproduce the 
old Dalitz result 



vf OA ---In^- - + 0(1/^) =0.01185, 
which should be compared with the present experimental value PP: 

pexp^^O ^ e+e-7) 



r(^o ^ 27) 



0.01213 ± 0.00033. (4.5) 



As already mentioned in the Introduction, the traditional radiative corrections to the total 
decay rate are tiny. The corrections corresponding to the first and second term in the decom- 
position of the QED corrections ()4.2() were first numerically evaluated in [TTI and analytically 
in [I2I with the result: 



r(7rO ^ 27) Vvr/ V9 m 9 m ^ 81 27 

= 1.04 X 10"''. (4.6) 
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Let us note that this formula is not based on the soft photon approximation, but the whole 
energy spectrum of the bremsstrahlung photon is included. A real photon emitted from 
the pion vertex is not considered. If we take into account the remaining NLO I'jR (ChPT 
corrections) and IjIR electromagnetic corrections we get an additional contribution 



-nNLO 
ChPT+l7/i? 



e ' e 



-7) 



27) 



(1.8 ±0.6) X 10" 



(4.7) 



where the error stems from uncertainty of C2-^^ Clearly these two corrections are small in 
comparison with the present experimental uncertainty ()4.5() . 

Similarly we could evaluate the corrected rate in the soft-photon approximation. The 
result then depends on AE. For AE ~ lOMeV, for instance, the result reads ° ~ 4x 10~^. 



4.4 Slope parameter 

We have now all the elements at hand in order to discuss both the extraction of the slope 
parameter from the data, and the prediction that can be made for it in the framework of 
the low energy theory. With the help of the formulae (|3.1U|) , 1)0.21^ , (|D.2|) and the definition 
(|3.6|1 we easily find 



M2 



967r2F2 



1 + ln 



1 /a 
360 W 



(4. 



where the individual terms in the square bracket correspond to the counterterm, the charged 
pion chiral loops, and the charged pion vacuum polarization function contribution, respec- 
tively (the latter we include here only for the sake of completeness, numerically it is negligible, 
being of the order 10~^, and thus can be safely neglected). Using the previous inputs, we 
obtain the following theoretical prediction for the slope parameter 



0.029 ±0.005. 



(4.9) 



As we have noted in the preceding section, previous experimental analyses, as a rule, did 
not include the contribution of the two-photon exchange (which was treated as negligible 
due to the superficial arguments based on the Low theorem). Therefore, according to the 
formula (|3.6j) . the systematic bias due to this omission can be roughly estimated as, cf. 



2 IhT 



- (- ) (27r2 - 3) = 0.005. 



8 Vvr 



(4.10) 



This corresponds to a shift of the central values for extracted from the Dalitz decay 
measurements which goes into the right direction towards the independent CELLO result. 



5 Summary and conclusions 

The present work provides a detailed analysis of next-to-leading order radiative corrections 
to the Dalitz decay amplitude. This study involves the off-shell pion-photon transition form 

^"Note that the ratio ^%^^1^2-~i)'^ independent of the unknown constants Ci and K.F to the order consid- 
ered here. 
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factor, which requires a treatment of non perturbative strong interaction effects. We have re- 
hed on representations of this form factor involving zero- width vector resonances. In contrast 
to the simplest vector meson dominance representation, our approach satisfies various short 
distance constraints from QCD. Our analysis also includes the one-photon irreducible contri- 
butions, which were usually neglected. We have shown that, although these contributions are 
negligible as far as the corrections to the total decay rate are concerned, they are however 
sizeable in regions of the Dalitz plot which are relevant for the determination of the slope pa- 
rameter of the pion-photon transition form factor. We have also obtained a prediction for 
which is in good agreement with the determinations obtained from the (model dependent) 
extrapolation of the CELLO and CLEO data. The present difference with the central values 
directly measured in the latest Dalitz decay experiments can be ascribed to the omission of 
the radiative corrections induced by the one-photon irreducible contributions. Unfortunately, 
the experimental error bars on the latest values of a-j^ extracted from the Dalitz decay are still 
too large to make a comparison with the CELLO and CLEO values meaningful. Nevertheless, 
we think that a precise measurement of a^r which would not rely on any kind of extrapolation 
remains an interesting issue. Hopefully, future experiments, like the one proposed by the 
PrimEx collaboration at TJNAF, will improve the situation in this respect. 
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A Form factor projectors 

Here we list the projectors A^ which allow to obtain the form factors F = P, A± and T from 
Eq. ()2.7|) . One possible choice is 



1 . „.,^fl, Tn 



^ [2A^-{{k-qf-{k-6f)5']Af^j,, (A.l) 



= ^^ie'^'"'^k,p+^p,pTY^[2A' + q\k-iqTdm.6m^ik.q)-q'^^]j, 

-^^^[4A2 - 2{k ■ {q T mk ■ q)S' +m^k-{qT S)nj^^k ■ 6 - 5^|^]75 
-^^[4A2 - 2ik ■ (q T S)){k . q)5' +m'{k-{q^ 5))']A^75, (A.2) 
= -^ie^'"'^k,p+^p.p, (A.3) 
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where we have introduced the shorthand notation g = + p_, A = (k ■ p+)p- — {k ■ p-)p-^., 
and 6 = p-^- — P-. In the hmit m ^ these expressions simphfy to the form 

= ^^■^e'^''''^hP+aP-p-^^[^'-2ik-p+){k.p_)6']A^^j,, (A.4) 
= TY^[2A^ + qHk.iqT6))ik-6)]h^ik.q)-qf^^]j, 

-^^^[4A2 - 2{k . {q T 6))ik ■ q)5'][rk ■ 5 - S'^^h^, (A.5) 

= ^ie^'^^^Kp+aP-p. (A.6) 



B The pion-photon-photon vertex 



As we have seen in the main text, the doubly off-sheh 7r^-7*-7* vertex, defined as 

I d''xe''-^{0\T{r{x)fm\7:\P)) = -ie^^^Pl^ppA^.^*^,{l\ {P - if), 

is a necessary ingredient for the calculation of the Dalitz decay amplitude. While in the 
case of the one-photon reducible contribution it is sufficient to use the corresponding form 
factor ^^0^*^* (0, i^) for P < M^, which is the region where ChPT (with virtual photons) is 
applicable, for the leading one-fermion reducible and one-particle irreducible contributions it 
is necessary to know A^o^*^*{P , {P — l)"^) as a function of the momentum I in the full range 
of the loop integration. In the present Appendix, we neglect temporarily the electromagnetic 
interaction, so that ^^o^.-y*(/^, {P — l)"^) will refer to the strong matrix element. We briefly 
summarize some basic properties of the form factor ^^o^*^. (/^, (P — l)"^) that are general 
consequences of QCD, as well as the results of [23) which we shall use in the sequel. The low 
energy expansion of ^^o^«^« in the presence of both strong and electromagnetic interactions 
will be the subject of the next Appendix. 

General properties of the form factor ^^o^.^. (/^, (P — l)^) in QCD were investigated earlier 
within various approaches j4()| 141 j. In the chiral limit, the on-shell value of the form factor is 
entirely fixed by the QCD chiral anomaly. Therefore, within ChPT the low energy behaviour 
is expected to be 



.A * ■ 



Nc 



l + O 



P-l t 



rrin 



AI^'A^^'Ah 



(B.l) 



where the higher order corrections come from pseudo-Goldstone boson loops, as well as from 
higher order contact terms. In particular. 



A 



(0,0) 



Nc 



l + O 



Ah 



Another exact result, the leading short distance asymptotics for / — 
from the operator product expansion (see IIB])- For A ^ cxd we have 



(B.2) 

oo, (P fixed), follows 



iiXlf,iP-Xl)') 
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+ ... 



(B.3) 



26 



The ellipsis stands for higher order terms in the short-distance expansion, or for 0{as) QCD 
corrections to the terms that are shown. Notice that the latter are not affected by quark 
mass effects, so that Eq. ()ij.3p holds beyond the chiral limit. On the other hand, the ex- 
pression HB.3|) assumes isospin and CP invariance of the strong interactions. The explicit 



form of ^^o^*^*(/^, {P — /)^) in the intermediate energy range, however, is not known from 
first principles. Among the various approaches that have been considered, models inspired 
by the large-A'^c properties of QCD have been proven particularly useful in order to pro- 
vide parameterizations of the form factor At^o^*^*{P , {P — l)"^) compatible with the above 
low and high energy behaviours predicted by QCD. Let us give here a brief overview of the 
results obtained in Ref. |27j within this framework. At leading order in the 1/Nc expansion, 
) can be expressed as infinite sum of the tree- level exchanges of the zero- 
width resonances in the various channels. Truncating this infinite sum and keeping only the 
contribution of the lowest resonances, i.e. the lowest vector meson octet in the present case, we 
obtain the Lowest Meson Dominance approximation {LMD) to the large- A'^c expression |26) 

aLMD (j2 ( p ,x2n _ P + {P - If + Ky 

^^0^,^, [L ) - ^ _ M^){{p - ly - m2 ) • ^""-^^ 

This Ansatz satisfies all the properties of ^^0^.^* discussed so far, provided the constant Ky 
is chosen such as to provide compatibility with (|B.1|) . 

= Aorr (0' 0) = - + Ci + • • • )• (B.5) 

The second equality involves the leading quark mass corrections described by the combination 
of low energy constants given in Eqs. (jClSf) and HC.14() below, while the ellipsis stands for 
higher order quark mass corrections, that will not be considered here. Let us note that if the 
large- A(7 vector meson mass is identified with the physical mass of the p meson. My = Mp, 
the form factor A^o'^P , {P — /)^) contains Ci as the only free parameter, and interpolates 
smoothly between ()BT|) and HB.3|) . On the other hand, at low energy, the non-analytical 
contributions from Goldstone boson intermediate states are not taken into account (note that 
according to the large- Ac counting rules, meson loops are suppressed in the I/Nq expansion). 
As further discussed in [221, simple Ansatz (|B.4|) is not sufficient to describe the full 
asymptotic behaviour for — > 00, where = — (fj'i+Q'i) ^^^^ fixed u = {qf — / (Qi ~^ Q2) = 
±1, given by the general formula j41j 



-4.o,*,*(gf,gi) = -^^ + o(ij) , (B.6) 



with a function /(w) that is not known explicitly. In order to reconcile the large Nc ansatz 
with HB.6() . at least one additional vector resonance is unavoidable. We thus obtain, in the 
notation of j27j, the more general Ansatz 

almd+V(^2 9. ^ gig^gi + gj) + + <i? + ^2^192 + ^5(gi + gj) + ^7 7^ 

^.0,.,. Wi,g2J 3 (,2_M2j(,2_M4)(,i-M2^)(gi-M2.) • ^^-'^ 
The chiral anomaly fixes now 

nr = A'Vr (0' 0) = - ^^Tf^^^ + Ci + • • • ), (B.8) 
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while the large behaviour of ^^o^*^* (Q^, 0) requires k\ = 0. From experimental data one 
can also determine 

K5 = 6.93 ±0.26 GeV^ (B.9) 

(one takes My^ = 769 MeV, My^ = 1465 MeV, = 92.4 MeV, further details can be found 
in ) . Finally, as pointed out in Ref . , the coefficient K2 is also available from Ref . |12 , 

K2 ~ -4(M^^ + Afv-J = -10GeV^ (B.IO) 
a value which lies within the range considered in Ref. j27j . 



C Chiral Expansion of the pion-photon-photon vertex 

In this Appendix, we first summarize the results of our recalculation of the pure QCD form 
factor At^o^*^* (0, in two- flavour Chiral Perturbation Theory [^[20113100] up to one loop, 
i.e. up to the order 0{p^). After that, we describe the additional modifications that appear 
if electromagnetic effects are also included. 

C.l e = 

The relevant chiral Lagrangian can be written in the form 

where the terms with even intrinsic parity at order O(p^) and O(p^) are 

£(2) = ^(^D^^u+D^U + x^U + U+x) (C.l) 

+^(x+t/ + U+x? + ^jiD.UD'^x^ + D^U+D,x) 

-h-(x^u-u^x? 

16 

+\{hi + h3){x^x) + ^(^1 - /i3)(detx+ + detx) 

+ h{R^,-L,,){R^'^-Ln). (C.2) 

The odd intrinsic parity Wess-Zumino-Witten Lagrangian, which accounts for the two-flavour 
anomaly can be written in the form (44j 

(C.3) 
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In the above formulae, the notation is as follows: 



( 7r° \/27r' 



D^\J = d^U -\r^U + \Ul^, ^^ = U+d^U, (C.5) 

Rfiiy = df^ri, - dyr^ - i[r^, r^], L^^, = (9^/^ - 9,,?^ - i[/^, (C.6) 

1 - 1 

= r^--(r^), If, = - ^{Iti), (C.8) 



rf, + l^i), v^y = df,Vu - dpV^ - i[v^,v„]. (C.9) 



Further relevant chiral invariant Lagrangians of order 0{p^) and also the other details can be 
found in [451 146j and references therein. 

The form factor ^^o^.^. (0, starts at the order 0{p^) with the tree graph with vertex 
derived from the Wess-Zumino-Witten Lagrangian ()C.3|) . and reproduces the anomaly result 

(ED), 

-4^o%(0,/^) = -^, (C.IO) 

since Fq can be identified with at this order. At the next-to-leading order, there are two 
types of one- loop contributions with one vertex from C^^^^, namely the tadpole and the 
bubble graphs (see Fig. EB- Another type of contributions correspond to the contact terms 



I I 
I I 



Figure 10: Next-to-leading order corrections to ttj'j vertex. 

derived from the tree graphs with one vertex from the odd intrinsic parity part of C^^^ . A last 
contribution comes from the renormalization factor of the external pion leg; this one contains 
the tadpole with vertex from and contact terms with vertices from C^^\ see Fig. 1111 
Putting all these parts together we obtain the following result 

0-- = + + A 

Figure 11: vr'^ wave function renormalization in next-to-leading order. 

A^o,.,,{0,f) = -Y^{l +Ci - ^[(i^)' (5 - 16vrV^(/^)J.(/^)) -C,]}, (C.ll) 



where crp(s) = Wl — AMp/s and 



Ms) = -1, r (c^i2) 
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is the Chew-Mandelstam function (the scalar bubble subtracted at s = 0). In the above 
formula, we keep the neutral and charged pion masses equal (within pure QCD their difference 
is an effect of second order in the isospin breaking parameter (m„ — m^)), the isospin breaking 
QED corrections, which are, taking e = 0{p), of the same order as the leading order terms, 
will be taken into account in the next section where further details can be found. Finally, Ci 
and C2 represent the following renormalization-scale independent combinations of the 0{p^) 
low energy constants Ai and Cj introduced in [35] and [3^] respectively: 



vr 



{A2 — 2A'^ 



64 



-vr 



C2 



3 



-11 



4cS 



on 

AAl)M^ + -^iAl + 2Al)2mB 



Ac'j)Ml + -{hcl + + 2cl)2mB 



In- 



Ml 
^2 



with m = (nid — my) /2 and 



Clip) = -QA^^MliA^ 



2^r 



13- 



(C.13) 



(C.14) 



The renormalization of the external pion line is responsible for the replacement of the constant 
Fq with the physical decay constant i<V in the leading order term of ()C.11|) . 

The actual values of the constants Ci and C2 are not known from first principles. Recently 
the relevant combinations of the low energy constants that occur in C2 have been estimated in 
|27j by using the matching of the LMD approximation to the large Nc form factor Almd{^''^ 
and the large Nc approximation to the ChPT result A'^'^^^{P,0). Since in the large- A'^c limit 
the contribution of meson loops is suppressed, the chiral logarithms as well as the running of 
the renormalized couplings with the renormalization scale are next-to-leading order effects. 
Following 1321) we assume that the values of the low energy constants obtained this way 
correspond to scale given by the mass scale of the non-Goldstone resonances n ~ My- We 
thus have the following LMD determination of the low energy constant C2 |27j : 



C5(Mi 



V)LMD 
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1 /47rF^^2 
AN^y Mv 



(C.15) 



KMyJ 

The same procedure can be done with the LMD + V approximation; in this case we find |2 



ANc\Mv, ) \MvJ M^M^^ 



(C.16) 



The issue of the quark mass corrections to A^^o^*^* (0, t^), contained in Ci, have been addressed 
in |3H]) and more recently in (22] • From (HHI, one infers 



Ci 



rud - rriu 



— m 



(0.93 ±0.12) ±0.14- 10 



-2 



(C.17) 



with fh = {niu + rrid) /2. 
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Numerically, with {nid - mu)/{ms - m) = 1/43, Af^ = 135 MeV, = 92.4 MeV, 
My = My, = 770 MeV, My^ = 1465 MeV and K5 given by (lR9l) . we then have the following 
determinations 

Ci = (2.2±0.3) X 10-^ 
CI{Mv)lmd = (1.5±0.5) X 10-\ 
CI{Mv)lmd+v = (1.8 ± 0.6) X 10-1. 

A 30% uncertainty, typical for a result based on a leading order large- A'^c calculation, has been 
assigned to C2. Within these error bars, the LMD result is stable with respect to the inclusion 
of a second resonance. Notice also that a variation of the scale between, say. My = My, and 
My = Myj gives C5(My,) — C2(My2) = 0.02, which is well within these error bars. 

C.2 e^O 

In this section we shall describe the results of our calculation of the next-to- leading 0{p^) 
corrections to the leading order amplitude in the expansion scheme in which the electric 
charge, fermion masses, and fermion bilinears are assumed to be counted as quantities of 
order p. Within this scheme, the O(p^) Lagrangian reads 

£(2) = _^(^Df'U+D^U + x^U + U+x)+e^ZF^{QUQU+) 
--^F^uF^"" + V^(i7 • {d - ieA) - me)iP, 

where 

Q = diag(2/3,-l/3) 

is the quark charge matrix. Then the leading order amplitude, which corresponds to the tree 
graph with one vertex from ^^^)zw^ ^® of the order 0{p^). Let us also note that electromagnetic 
splitting of the charged and neutral pion masses is treated as a leading order effect. The 0(p'*) 
Lagrangian with even intrinsic parity then reads 

/'(4) _ r (4) I £(4) ,£(4) 

The C^p^ is the same as ()(].2|) while the explicit form of >Cgtp2 and ^\tpt found in [20] 

and In the following we need only 

{kl{D^'U+D,,U){Q^) + k2{D''U+ D,,U){QUQU+) 
+h {{D^'U+QU){D^U+QU) + {D^'UQU+){D^UQU+)) 
+k4{D''U+QU){D^UQU+) + ...}, 

The NLO contributions within the pure QCD were presented in the previous section. In 
the enlarged case there are two main distinctions: First, because the pion mass difference is 
of order O(p^) now, we have to take care of the unequal pion masses in the loops. Second, at 
order O(p^), the pion decay constants F^o and F^± are different as a consequence of the new 
type of 0{p'^) terms coming from 21 as well as of unequal tadpole contributions. The mass 
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difference leads to additional terms of the form In (M^±/M^o)) the latter to the replacement 
of the Fq with F^o (and not with F^± ) in the leading term as a result of the renormalization 
of the external pion line. Taking all these effects into account leads to 
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(i + ln^-167rV2+(/2)j^+(/2)^ (C.18) 



with Ci and C2 now given by 
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Because the constant -Fj^o is not known very accurately, we use the following relation [241 125j 



F^o=F^(l 
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with 



}Cp = {3 + ^Z)ki - fZk2 - 3k3 - AZki 



(C.19) 



(C.20) 



to eliminate F^o in favour of F^, where F^^ = F^± |e=o is measured in the charged pion decays, 
see Pni- We thus write 



A 



ChPT (i2 



where 



647r2 



)Cf + 



Nc 



127r2K, 



In- 



Wir'^F^ M 



M2 



167^2 F2 (3 



+ In 



M2 
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The fej, i = 1, ... ,4 are the a priori unknown scale-independent constants, defined in terms 
of the bare low-energy constants from C^^^2 according to the formulae 
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where 

27 1 3 
^1 = - ^2 = = 2Z, as = --. (C.22) 

In order to obtain a numerical evaluation of /Cj? we further express it in terms of the 
analogous SU{3) constants Ki. Several determinations of the latter are available in the 
literature j37j . In the most recent work I5H|, a new estimate of these parameters based on 
sum rules involving QCD 4-point correlators (for SU{3) case, parametrized with help of the 
improved chiral Lagrangian with resonances in the spirit of large-A'^t; approximation) was 
made. In order to use these results, we first have to match the SU{3) variant of the theory 
with the SU{2) one we have used in our calculation. This can be made as follows. Starting 
from the fact that ICp enters the formula HC19|) expressing the electromagnetic difference 
between F-^ = F^±|e=o and F^o, we can write (in the SU{2) case) 

!^F = l-^--r^,ln-^. (C.23) 



647r2 ^ 167r2F2 M\ 

The ratio F^a/F-,^ can be calculated within the SU{3) version with the result^^ EH 
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F^ 



= l-^ln^-^(41n(4)+lnf4l)+l 
SU(Z) 167r2F2 327r2 V \ ) \ fi^ J 

+ {^K[ + ^K^, - 2KI + Kl + ^^Kl + 



therefore, upon matching the two expressions, it follows that 
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jKf = 1 



647r2 ^ F^ 
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5C/{3) 167r2F2 M2o 



ln-^= (C.24) 



Inserting into this expression the values Kl = —2.71 x 10~^, K2 = 0.69 x 10^'^, = 
2.71 X 10-3, Kl = 1.38 X 10-3, Kl = 11.59 x 10-^ and K^ = 2.77 x IQ-^ at a scale = 770 
MeV, obtained in PS from the lowest meson dominance approximation to the large- A^c" limit 
of appropriate QCD correlators we find 

/Cf = -28±8. 

Again, we have assigned to this value an uncertainty of 30 %, typical for calculations based on 
the leading order in the large- A'^c expansion. Although ICp '^s scale independent, the estimates 
of the low energy constants K^'lp) it involves depend on the scale at which they are identified 
with the resonance approximation. Varying again this scale between the values My = Mvi 
and My = My-j induces a variation in Kp which corresponds to these same error bars. 



^^The known experimental value of F^o = 92 ± 4 MeV has an uncertainty too large to provide a useful 
determination of ICf- 

^^Let us note, that the 5(7(3) on-shell 7177 amplitude .iO contains (besides the electromagnetic difference 
between and F^o) an additional 0{e^) contribution which originates in the electromagnetic correction to 
the Tvi] mixing. Within the SU{2) power counting it is in fact of the order 0{e^p^), so that it need not to be 
included in the matching procedure. 
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D NLO corrections to Fi, F2, H 

D.l Corrections to the vacuum polarization function 

The vacuum polarization function n(Z^) starts at 0{p^) with three types of contributions, 

n(z2) = n,±(/2) + n,±(/2) + ncr(/'), 

where the first two correspond to the pion bubble and tadpole, and to the fermion bubble, 
respectively, and the third one is a contact term from jC^^\ which is necessary to renormalize 
the UV divergences. In dimensional regularization one has 



n..(.) = n..(.)-^(in^ + i)-£ 
n,±(.) = n,±(s)-^(in^ + i)- 



^-i(ln4.-7+l) 



2a 
Stt 



1 



d-4 



{In An -7 + 1) 



(D.l) 



where 11 are the corresponding quantities in the on-shcll rcnormalization scheme in which the 
finite part of the counterterms is unambiguously fixed by the condition 11(0) = 0. We have 
then the standard result 



TT 12 74. 



dx<T^±(a;) _ a 



4^2^ X X - s 



IStt 



l + 247rV^±(s)J,±(5) 



(D.2) 



and 



n,±(s) 



a 



as f'°° dx ae±{x) / 2m'^\ 

TT 3 74^2 X X — S \ X J QtT 



1 + AStt^ 1 + 

V s 



(D.3) 



In our notation the 0{p^) counterterms contribute as 

ncr(s) = 16Tra(2h2 - - xs^ , 



where 



h2 
hi 



L__l(ln47r-7 + l)). 



2 1 



3(47r)2 \d 



While h2 renormalizes the divergent part of n^±(s), xg does the same with Tl^±{s). Note that 
in this scheme 11(0) 7^ and, as a consequence, rcnormalization of the external photon line 
has to be included by means of the factor 



z]/^ = 1 - in(o) 



1 - STra 2hlii,) - Mli) - xli,,) + _ In + 1 + — In 



67r 



m 



a 



2Att 



+ 1 
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D.2 Corrections to the fermion self energy 

In the same way, we can recalculate the fermion self-energy 

E(g) = ^Sv^(g2) + S5((?2), 
with the loop and counterterm contributions 

where (to regularize infrared divergences we have to introduce virtual photon mass m^) 



1 xim? - w?) - x{l - x)q^ + a { ^ rr? 3 

Q \1 ) — -in— I dxln 5 — m— \\n.—^ + - 



SrV) = £/ dx(l-x)ln - - - ■ ■■^ +--(^ln^ + 2 



and 



= -47ramx7, = - ^^(^3^ - ^(ln47r - 7 + 1)), 

E^((?2) = _47rax6, = + - ^(ln47r - 7 + 1)) . 

From these formulae, the fermion mass renormalization follows 

m = me + S5(m^) + mS\/(m^) 

= me -m |^47ra(4(/x) + + ^ (^^^^ ~ ' 

where m is the physical fermion mass. For the fermion wave function renormalization we 
need 



-|^=r„ = 2mS^(m') + Sv(m') + 2m'sV(m' 

a /I , vn? 1 m^ . t ». . ^ 3 
^ + - In ^ 

i2 ,,2 



TT V2 m2 4 /x2 ''^'^^ 



Thus, one has 

r^-i -, dT^iq) . ^ o: / 1 . vr? 1, m^ , 1 r, \ 3\ 
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D.3 Corrections to the form factors Fi_2 

We have the following standard formula for M^o > s > 4m?, 

^2 r /l_ae(s) 



a m 



In 



+ ivr 



(D.4) 



and 



^ la/1 1/, , af3 1, m? „— , ^ 

.2 5 



167rV, 



■is))} 



7rV2 AV2 m?, / mH 



+ - 



1 



1 



2 7o ~ sx{l — x) 
The counterterm contribution is 



In 



— sx{l — x) ■ 



We can compare now 

= F.(0) = 1 + - i in - i in 5 + 4.^..J(,)} = 



a r 3 1 . 1 . 

where the last identity is a consequence of the Ward identity. 



D.4 Complete l-yR LO+NLO form factors 

Putting the results of the previous subsections together one obtains (note that we must include 

1/2 

the external line renormalization factor Z-y Z^) 



^ljR;LO^^^ y) + A^^^'^^^^X, y) 



Fi{xMl,)-Ii{xMl,) 



+ (ZV2 _ 1) + (z^ _ 1) + a%\''J{xMlo) 
e^Nc i 



127r2F^o xM\ 



2m{Fi{xMl,)-Il{xMlo) (D.5) 



- 1) + {Z^ - 1) + -^F^ixMl,) + a^fJixMl,)) 
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Let us write U{s) = U{s) - n(0) = n(s) + {Zy - 1) and define = 1 + Fi(s) - Fi(0) 

Fi{s) + {Z^ - 1). Explicitly, for M^^ > s> Am^, 



n(s) = n,±(s) + n,±(s) = -{i + 



2 - , 1 



+ 



1 



1 - (Tejs) 

3scJe(s) "'Vl + aeis) 



■In 



Fi(s) = l + -|-1 + 



SaeiS L 



s — 2m 



-s) In 



s 6 
8m^ 



(T^±(s)|cT^±(s)| arctan 



1 



+ 



ivr 



3sc7r{s) 



s — 2m 



8m^ 



(D.6) 



1 - CTeis) 
l + fTe(s) 



(s - 2m^ 



1 /l-'Te(s) 



In 



In 



1 - q"e(s) 

l + fTe(s) 



+Li2 



0"e(.s) - 1 
2cTe(s) 



1/1-0^ 
2 V 2(Te(s) 



In 



1 + cre{s )\ it' 



+ 



ivr 



,2™^-2.)-i(,-2»^)l.,(i^ 



SO-e(s) - 

+ — In^ 1+ s-2m2 

27r \m£/ I ^scTpl 



In 



)]} 

1 - cre{s)\ ^ .J I 



+ ITT 



(D.7) 



Then, taking F2{s) = F2{s) and introducing a physical charge e = eZ}/^"^ (where e^/(47r) 
a = 1/137, . . .), we can rewrite (|D.5|1 in the form 



T^^^'^{x,y)+T^^^'^^{x,y) 



e-^Nc 1 i ^ 



e^Nc i 



-F2(xM2o), 



Utt^F^ xM% I 



2m(Fi(a;M2o) -n(xM2o^ 



Identifying now the leading order amplitude with the substitution Fi = 1, F2 = IT = a^^j^ 
in the above expressions, and using (|3.7j) we obtain 



rljR 



+ 



2xMl, 



and 



^NLO 



(x) = 2Re^Fi(xM^o) - n(xM^o) + a%fJ{xMl,) + ^ M\x + 2m^ ^^^''^''-'^ ~ ^ 



M2ox(1 + y2) + 4^2 



F^ixMl,) - 1 



Before concluding this section, let us give a brief survey of the IR divergent contributions. 
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They can be extracted from the formulae given above and read 

,2X ( 1 



a 



Fi{s)iRdiv = — In ^ W 1 + (s - 2m 



27r 



m 



m. 



, , 1 — Geis) 
In ( 1 + iTT 



F2{s)lRdiv = 'n.{s)lRdiv = 0. 

Thus, the IR divergent parts of the form factors are 
SPiRdiv{x,y) = 0, 
SAjjidivix,y) ■ 



1 — 



1 f ^'^g 



Fi{xM^o)iRdiv 



Inserting these expressions into formula (|3.7|) yields, after some simple algebra, 



-'nlo 



IRdiv 



m 



mi 



In ^ 1 + 1 



2rn^\ 1 / l-ae{xMX 



(D. 



E Loop functions 

This appendix is devoted to the so-called Passarino-Veltman [3^] one-loop integrals used in 
the main text. Generally one defines (working in d dimensions) 



i-K'^To{n) = (27r/i)^ 



(E.l) 



It is, then, common to denote these n-point functions in alphabetical order, i.e. instead of 
T one uses for 1-point integral the symbol A, for n = 2 ~ the B and so on. For the scalar 
functions and special combinations of arguments needed in our work we get successively 



-^-i(ln4vr-7 + l) 



2 

m 
In — ^ 



1, 



(E.2) 



^ ^(ln47r-7 + l) 



d-4 2 

^ -(ln47r - 7 -I- 1) 



d-4 2 



In ^ + 167r^ Jom(m^) 

2 2 

, m ( m 



In 1 



2 



(E.3) 



Co (0, m\. , ; , , 0) 



vr 



2 6Li2(^ + ie) 



6(m^ — m?) 



(E.4) 



Notice that according to this definition the loop functions are renormahzation scale dependent and conse- 
quently the bare LECs are also scale dependent. 
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Re Co(m^ M^, m^; 0, 0) 



+ Li2(l 



Li2 1 



2VA 



— — m? 



+ Li2 {\ + 

-Li2(l- 



2\/Am2 



{m\ — m2)('\/A + — m?) — {m\. + m?)M'^ 

(m^ — m?){-\/X + m\ — m?) — {m\ + w?)M^ 
2\/Am^ 

{m\ — m?){^f\ — m\ + m?) + (m^ + m?)M'^ 



(E.5) 



-I 

6 J 



with A = \{Ml,m\,w?) = (M^ — m? — m^)^ — Am^vn^, and m\ = m? + 5m^, where 
The four-point function appearing in Section 4.3 is given by: 



ReDQim^, 0, m^, iW^, mi, m?.; 0, m^, m^, 0) 



2y 



M'^m?{y'^ - 1) 



X <^ lo£ 



(ml — m?)(im^_ — m^ 



lo, 



igy + Li2(l-y)-Li2(l-y-i)}, (E.6) 



where y = ^{—b + Vb"^ — 4ac), with 



a = c = 



m^ 



m 



■ ((m+ - m^) (m?. - m^) + 2M^m^) . 



Asymptotics of the loop functions for A; ^ (x ^ 1), m fixed, read: 



Bo(m^;0,m^) = 2 + Bo(0,m^,m^) - 



5mi L /Smi 



In 



m^ 



+ ivr 



O 



5ml \ 2 



m^ 



Co(0, m^, m^; m^, m^, 0) 



In 



V m^ / 



+ i7r- 1 



1 5m,t 



4 m^ 



/5ml\2^ 



Co(mlM2,m|,m2,0,0) =Co(mlM2,m2,mlO,0)--^^ In(^) +0(1) 

Mi L V m^ / 



Do{m\0,m\Mlml,ml,0,m\m\0) = -1^ ln(^^) + 0(1) 
Asymptotics of the loop functions for m ^ 0, 5m1 > fixed, read: 



ReCo(0, m^jm^-^m^ jUi^ , 0) 



Re So(m2^; 0, m^) = -2(^-i(ln47r- 7+ 1)) +l-ln(^) +0(m2), 

1 



Smj. 
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ReCo(m^ M^, m|, m^ 0, 0) 



M2 - 6ml 



2Li2 1 



1 , 

+ - In 
2 



M2 



2 In 



2 

m \ 
Mi) ~ 



In 



± 



ReDo{rn'^,0,m'^,M^,m\,m^_,0,m^,m^,0) 



+ 0(m2) 



(5m2 L2 Vm2; 



1 /m^ \ f5m?,6m?_ 



+ - In^ 



vr 



-+0(m2) 



F Soft photon singularities 

In this Appendix, we briefly address the question of soft photon singularities, which are of 
relevance for the discussion in Section 3.3. We wish in particular to elaborate in somewhat 
greater detail on the statement made at the beginning of Section 3.3, concerning the absence 
of contributions that are independent of in the difference 7W^o^e+e-7 ~ 
the present context, we may arrive at this result as follows. First, note that the Ward identity 
(|2.22|) can be solved by the expression^^ 

(where ae = ^2(0) is the anomalous magnetic moment of the fermion), which includes the 
leading and next-to-leading order singularities for /c — > 0. Indeed, for the combination 

e^'{k)*uAf,{p^,p^ + k)S{p- + k)v, 

it is not difficult to prove that, for k such that (p_ ■ k) ^ 0, with p'^ = m? and p- fixed, 

e''ikyuA^{p-,P- + k)S{p- + k) 

= e^iky ^^-^ - - (^-^ ^ - : - "^-^^^^ + + on) + ak, (p_ . ky. 

2(p_ • k) 

Here (and in what follows), the remaining 0{1) terms, which are not written explicitly, are 
independent of k. In the same way, for k, (p+ ■ k) ^ 0, p'^ = m? and p+ fixed, we find 



S{-p+ - k))A^{-p+ - k, -p+)ve^{ky 

2{p+ • k) 



ve^iky + Oil) + Oik,{p+-k)). 
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Of course, the minimal solution can be written in the form 



^ r^Oe-e+(P- +k,P+) - , r^Oe-e+(P-.P+ + fc)] 



which takes into account only the leading order singularity for fc ^ 
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Notice that if in the above expressions one restricts the vertex function A^{qi, ^2) to its longi- 
tudinal part given by Eq. ()2.17|) . one arrives at the same expression, but with Oe replaced by 
— Ey(m^), which is both gauge dependent and infrared divergent. Including the contribution 
from the transverse part A^(g',g) cures both problems, and yields the anomalous magnetic 
moment Og. This can be checked explicitly at the one loop level with the expressions available 
in Refs. [H] and [221 • 

Let us recall that the one-particle irreducible (semi-)off-shell ■K^-e^-e'^ vertices r^Oe-e+ {P- + 
k,p^) and r^Oe-g+ + k) are free of poles for {p± ■ k) ^ 0. The same is true for the 

one-particle irreducible (semi-)off-shell e"'"-e~-7 vertices A^(p_,p_ + fc) and A^{—p^ — k, —p+). 
We have therefore, for {p± ■ k), A; — > and p± fixed, according to (|2.25|) . 



2(p_ • k) 

+ h vP^o^-Mm ,m ) 

+0{l) + Oik,{p+-k)), 

where, as above, the implicit 0(1) terms are independent of k. From this formula we can read 



off the associated Low amplitude given in Eq. (|3.22|) which, according to Low's theorem 
corresponds to the leading singular terms in the expansion of the complete amplitude in k, 
ip± ■ k) ^ , (with p± fixed) in the sense that the /c-independent 0{1) terms coming from 
r^'^^ are in fact cancelled in the complete amplitude by the corresponding 0(1) 

terms from the Tjf^ (let us recall that the one-photon reducible amplitude is of order 0{k)). 
On the other hand, we have^^ 

A^Sri-^ = e''{krur]f^'P'''% = M^SZe^^^, + + 0((p+ • A:), k). 
Therefore, the following subtracted quantity 

plV;i?,res ^ pl^^iJ _ ^HR,pole ^ ^^-^^ ^ ^^^^^ . 

is both transverse and with the 0(1) terms independent of k. It can thus be expressed in 
terms of form factors P, A±, T, see (|2.f)|) . 

r^'^^'-f (p+,p_, k) = P^^^^'^^ix, y)[{k ■ p+)p^_ - {k ■ p_)p(t]75 

+ Al^^^'^'^x,y)[^p^^ - {k-p^)rH-A'_^''''-'^x,y)[^p^_ - {k ■ p^h^^]^, 
-iT'^'''^-a{x,y)a^''k,j,. 

Because these form factors are in fact 0(1), i.e. y]^^'^'^^ = 0{{p+-k), k) and also t]2^ = 0{k), 
we may conclude that the contribution of ^j^^'^^^ is tiny (it is suppressed by a factor a w.r.t. 
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Here we use the identities 

- 7m + -z—aecr,i,^k 



2{p- -k) V 2m 

and 



S-'(p_ + fc) = (1 - Ev(m2))((|>_ +t)-m) + 0{p- ■ k). 
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r^''^^) in the full kinematical region z/^ < x < 1. On the other hand, we should expect that 
the remaining gauge invariant combination, namely 

Mi^ie^,-, + m^j:::.^ = M^jZe^,.^ + oik, (p± • k)), (f.i) 

might be important for x sufficiently close to one (i.e. A; ^ 0) in spite of the suppression by 
a factor a. In the formula ()F.1|) . the one-particle irreducible part of the amplitude 



Ml^,^, = urlf'{p+,p-,k)ve>^{kr 

corresponds to the photon emission from internal lines, being therefore of the order 0(1) for 
A; — > 0. Notice also, that the Low amplitude is transverse, therefore we can decompose it in 
terms of P^""", A^""" and T^"*" form factors, where 



{p^-k){p+-k) M4„(l-x)2(l-y2)' 
± ""m 2{k-p±) ^ mM^,il-x)il±y)' 



2{p_ ■ k) 2{p+ ■ k) 
^^M2„(l-x)(l-y2)- ^'-'^ 
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